Slow relaxation in granular compaction 
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Experimental studies show that the density of a vibrated granular material evolves from a low 
density initial state into a higher density final steady state. The relaxation towards the final density 
value follows an inverse logarithmic law. We propose a simple stochastic adsorption-desorption 
process which captures the essential mechanism underlying this remarkably slow relaxation. As the 
system approaches its final state, a growing number of beads have to be rearranged to enable a local 
density increase. In one dimension, this number grows as N — p/(l — p), and the density increase 
rate is drastically reduced by a factor e~ N . Consequently, a logarithmically slow approach to the 
final state is found p^ — pit) = 1/lnt. 



PACS numbers: 05.40. +j, 81.20.Ev, 82.65.My 



Systems consisting of many macroscopic particles such 
as sand and powder exhibit complex behavior despite 
their apparent simplicity jl| . Shaking sand may result in 
size segregation, richpattern formation m , or creation 
of convection rolls Adding a few grains to a sand 
pile leads to avalanches This wealth of phenom- 

ena has stimulated analogies to other intriguing physical 
problems including the motion of flux lines in supercon- 
ductors, metastability in glassy systems, and surface in- 
stabilities in fluids Q. 

Despite growing interest, the physical principles under- 
lying sand remain largely unknown. Although individual 
grains are solid, it is inappropriate to classify their collec- 
tive properties as entirely solid-like or liquid-like. Con- 
ventional thermodynamic theory is not applicable to sand 
as thermal fluctuations are negligible, fc^T = 0. Most of 
the experimental and the theoretical studies have focused 
on spatial inhomogeneities of one sort or another, and 
much less attention has been given to equally fundamen- 
tal but simpler situations such as compaction. Granular 
compaction can be viewed as a model system for non- 
thermal relaxation in a disordered medium. 

Granular compaction is relevant to production, pack- 
ing, and transportation of a wide array of industrial and 
agricultural products such as food, grains, chemicals, 
drugs, and building materials. A granular assembly pro- 
vides us with a practically uniform system where upon 
vibration, the well-defined bulk density evolves from a 
loosely packed mechanically stable initial state into a 
denser final state (see ref || and references therein). The 
system exhibits no convection, spatial patterns, or oscil- 
latory behavior. The system slowly explores available 
microscopic configurations, and eventually low-density 
metastable configurations are eliminated. Several mecha- 
nisms have been proposed to explain the kinetics of com- 
paction. The concept of compactivity, which quantifies 
how far the density is from its maximal value, was argued 
to play the same role as temperature does in conven- 
tional statistical mechanics [^0[-p!2[. Such a theory pre- 
dicts fast exponential relaxation of the density with two 
relaxation times associated with collective and individual 
modes [O. According to another theory, the motion of 



the voids filling the space between the particles is effec- 
tively diffusive, and as a result a power-law relaxation 
is predicted 0,^5) ■ Although the proposed mechanisms 
are compelling, their quantitative predictions fail to de- 
scribe the time dependence observed experimentally B . 

In the compaction experiment of Knight et al [p|, 
monodisperse 2mm diameter glass beads are confined to 
a 1.88cm diameter, lm long pyrex tube. The tube is 
tapped vertically by a vibration exciter. The duration 
of each pulse is much smaller than the waiting time be- 
tween different taps so that all of the kinetic energy of 
the beads is dissipated, and the beads are at rest before 
the next tap. The column density was measured by a 
noninvasive capacitive technique. The beads are initially 
loosely packed with a volume fraction of 0.577 ± 0.005. 
Precautions were taken to prevent humidity (the tube is 
evacuated), convection, and electrostatic charging. The 
vibration intensity is parameterized by T, the ratio of 
the peak acceleration of a tap to g — 9.81 m/s 2 . In 
general, the steady-state density increased monotonically 
with increasing T. In Fig 1, the time dependence of the 
packing fraction, pit), is shown for a representative vi- 
bration intensity S. Shown also is a four-parameter fit 
p{t) = pj — A/9oo/[l + i31n(l + </T)]. The parameters pf, 
Apoo, B, and r depend only on the acceleration T. This 
inverse logarithmic form fits the experimental data bet- 
ter than exponential or algebraic laws suggested by the 
aforementioned compaction theories. Logarithmic relax- 
ation has also been observed in the the decay of the slope 
of a vibrated sandpile ||. 

As the compaction progresses, individual grains move 
slowly, and when a grain size void is created due to such 
motion, it is quickly filled by a new grain. New grains can 
not move into space occupied by other grains. In other 
words, they interact with their neighbors via a hard core 
interaction. When the packing fraction is large, voids 
the size of a particle are rare and a large number of voids 
have to be rearranged to accommodate an additional par- 
ticle and a local density increase. Following this line of 
reasoning, we propose a useful heuristic argument. Let 
us denote the volume of a particle by V and the pore 
volume per particle by Vq. Thus the packing fraction 
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is p = V/{V + Vo), or alternatively, V = V(l - p)/p. 
Assuming that N particles are rearranged in such a way 
that they contribute their entire free volume to create 
a particle size void, NVo — V. We find N = p/(l — p). 
The time associated with such a rearrangement should in- 
crease exponentially with N, T ~ e N . Consequently, the 
density increases according to the following rate equa- 
tion dp/dt oc e~ N = e~ p " l ~ p K The solution of this 
equation is given asymptotically by p(t) = 1 — 1/ hit. 
We made the unrealistic assumption that all of the free 
volume is available for large void creation. Furthermore, 
we ignored the structure of the granular assembly, and 
the fact that the density cannot exceed the close-packed 
value. These features can be avoided by considering the 
one-dimensional situation where a void can be uniquely 
"assigned" to its neighboring particle. Interestingly, the 
exponential reduction factor agrees with analytic results 
in ID. 

0.620 | 




Packing fraction vs. time. Compaction data (o) near the 
bottom of the tube for acceleration V = 2.3. The data rep- 
resent an average over 5 different runs and the error bars 
correspond to the rms variations between runs. The solid line 
is a least square fit to the inverse logarithmic form discussed 
in the text. 

We now propose a simple adsorption-desorption the- 
ory for granular compaction. The corresponding stochas- 
tic process was previously studied in the context of 
chemisorption Jl6|-[l8| and protein binding jnj. The 
model is defined as follows: Identical spherical parti- 
cles of unit diameter adsorb uniformly from the bulk to 
a substrate with rate fc+ and desorb with rate k—. In 
other words, fc+ adsorption attempts are made per unit 
time per unit length, and similarly, the probability that 
an adsorbed particle desorbs in an infinitesimal time in- 
terval (t,t + dt) is k-dt. While the desorption process 
is unrestricted, the adsorption process is subject to ex- 
cluded volume constraints, i.e., particles can not adsorb 
on top of previously adsorbed particles. The attempted 
adsorption event in Fig. 2 is thus rejected. This stochas- 



tic process is well-defined in arbitrary dimension. How- 
ever, we restrict our attention to one-dimension where 
analytic results are available. This "car-parking" process 
clearly satisfies detailed balance, and after a sufficiently 
long time, the system reaches equilibrium. The proposed 
model can be viewed as a homogeneous one-dimensional 
gas of hard spheres [|o| . 

To analyze the parking model, it is useful to con- 
sider first the simpler lattice version of the process where 
particles occupy a single lattice site. In this case the 
density, p, satisfies the Langmuir mean-field equation 
dp/dt = fc+(l — p) — k-p. The gain term is proportional 
to the fraction of unoccupied space, while the loss term is 
proportional to the density itself. The steady-state den- 
sity, Poo, which is obtained by imposing dp/dt = 0, de- 
pends on the rate ratio k = k+/k- only, poo = k/(l + k). 
The system relaxes exponentially towards its final state, 
Poo — pit) — (Poo — Po)e - */ T , with po the initial density. 
The relaxation rate is simply a sum of the adsorption 
rate and the desorption rate t _1 = fc+ + 
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Fig 2 The adsorption-desorption process. 

On the continuum, however, only a fraction of the 
empty space is available for adsorption. It was previ- 
ously shown that in equilibrium the probability s(p) that 
an adsorption event is successful equals e~ p /( 1 ~' 3 ) [H. 
This so-called "sticking coefficient" equals unity in the 
low density limit, and the above mean-field treatment 
holds. On the other hand, in the high density limit this 
coefficient vanishes exponentially, s(p) oc e-W-P), as 
p — ► 1. The excluded volume interaction effectively re- 
duces the adsorption rate, k+ — > k+(p) = k+s(p). Thus, 
a modified Langmuir equation can be written 



This evolution equation was constructed to give the ex- 
act equilibrium density, and is reasonable as long as the 
system is sufficiently close to equilibrium. 

Using Eq. ([[]) , the equilibrium density is obtained from 
the following transcendental equation, ae a — k, with 
a = Poo/(l — Poo)- We find the following leading behav- 
ior in the two limiting cases 

Poo(fc) - {l-(lnfc)- 1 fc>l! (2) 
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While the mean-field linear relation is recovered in the 
dilute limit, the dense limit is characterized by a loga- 
rithmic cusp. By contrast, mean-field predicts a power- 
law dependence of the density p^ = 1 — 1/fc, for k 3> 1. 
The effect of the volume exclusion constraint is striking, 
a huge adsorption to desorption rate ratio, k = 10 9 , is 
necessary to achieve a 0.95 steady-state occupancy. It 
is also interesting to consider the equilibrium void dis- 
tribution. The density of voids of size x is exponential 
■Poo (a:) = (3e~ ax , with (3 = p^/(l - Poo) and the pre- 
viously defined a |Tq ] . Indeed, the sticking coefficient is 
proportional to the density of voids larger than a particle 
size, s(p) oc J 1 dxP oa (x). 

We now focus on the relaxation properties of the 
system. The granular compaction process corresponds 
to the high density limit, and we thus focus on the 
desorption-controlled case, k 3> 1. Hence, let us fix k+ 
and consider the limit k- — > of Eq. (|l|). The early 
time behavior is dominated by adsorption and can be 
obtained by neglecting the desorption term. For simplic- 
ity, we consider a vanishing initial concentration. For 
sufficiently early times, t <C the density grows lin- 

early in time, p(t) = k + t. At the later stages of the 
process, t I/fc+, the system approaches complete cov- 
erage, poo = 1, according to p8 21 



P(t) = Pc 



ln(fc+t) ' 



(3) 



This is confirmed by numerical simulations in one dimen- 
sion. Use of Eq. (|l|) is justified a posteriori since the sys- 
tem evolves slowly and has enough time to equilibrate. 
The inverse logarithmic behavior is simply a reflection 
of the exponentially suppressed adsorption in the dense 
limit. In writing Eq. ([?) we neglected higher order loga- 
rithmic corrections such as lnln(fc+i). We conclude that 
the excluded volume interaction gives rise to a slow re- 
laxation. 

Eq. (|^) holds indefinitely only for the truly irreversible 
limit of the parking process, i.e., for k — oo. For 
large but finite rate ratios, the final density is given by 
Eq. (||). As the density approaches this steady state 
value, the loss term becomes significant and should be 
taken into account. The crossover time between the two 
different relaxation regimes, to, can be conveniently esti- 
mated by equating the time dependent density of Eq. (|^) 
with the equilibrium density of Eq.(||) 1 — l/ln(fc+to) = 
1 — 1/ ln(fc+//c_ ) , and as a result to = l/fc_. For t 3> to, 
the loss term is no longer negligible. By computing how 
a small perturbation from the steady state decays with 
time, an exponential relaxation towards the steady state 



is found 
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for t ^> to- The relaxation 



time is indeed to, however, an additional logarithmic cor- 
rection occurs, r = t (l — Poo) 2 — to /(lnfc) 2 . The above 
results can be simply understood: the early time behav- 
ior of the system follows the irreversible limit of = 0. 
Once the system is sufficiently close to the steady-state, 
the density relaxes exponentially to its final value (Fig 3). 



Hence, two relaxation curves with sufficiently large final 
densities are indistinguishable over a significant temporal 
range. 




Fig 3 The coverage versus time. Shown are solutions of 
Eq. (1) with fc+ = 1. The upper curve corresponds to the ir- 
reversible limit poo = 1, while the lower two curves represent 
reversible parking with equilibrium densities of — 0.7 and 
0.8 respectively. 

In the desorption-controlled limit, the void left by a 
desorbed particle is quickly filled by an adsorbed par- 
ticle. Hence, the desorption mechanism can be seen as 
playing a similar role to diffusion of particles. Indeed, 
theoretical studies show that adsorption processes where 
adsorbed particles diffuse rather than desorb give rise to 
the same density relaxation as in Eq. p2| . Our treat- 
ment was restricted to one-dimensional processes, but we 
expect that the results hold in higher-dimension as well. 
Hence, the above relaxation is rather robust and insensi- 
tive to many microscopic details. 

The inverse logarithmic density relaxation towards the 
steady-state is the same one observed for granular com- 
paction (see Fig 1). Furthermore, the experimentally ob- 
served relaxation curves which correspond to large com- 
paction, r > 2.7, are indistinguishable over the observed 
time range consistent with the above theory. The largest 
packing fraction observed in the experiment are well be- 
low the maximal close-pack value in 3D. Thus, according 
to the parking theory the logarithmic relaxation regime 
may be followed by an exponential relaxation. There 
exists a maximal time scale which characterizes the tem- 
poral behavior. Such a prediction is consistent with the 
experimental observation that steady state is achieved 
within a finite time. Finally, once the system reaches 
its final state, density fluctuations can be measured and 
compared with Monte Carlo simulations of the parking 
process. Our preliminary results indicate that measured 
and simulated Power Spectrum Density (PSD) of the 
density fluctuations are strikingly similar p4j . The high 
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frequency limit is Lorentzian (/~ 2 ), while the low den- 
sity limit is white- noise (/ ). The intermediate regime is 
approximately a powerlaw f~ a , with 0.7 < a < 1.4. 

In a realistic granular material, an individual particle 
can not penetrate its neighbors, but it also must be in 
contact with several other particles. Our model properly 
accounts for the hard core repulsion, but it ignores me- 
chanical stability. We argue that in the long time limit 
mechanical stability can not play a significant role in de- 
termining the motion of individual grains during the com- 
paction process. Instead, the motion is limited primarily 
by the presence of other beads. The tradeoff, however, 
is that such a simple theory can not make predictions 
about the final density. Nevertheless, it does elucidate 
the leading mechanism in granular compaction. It will 
be very useful to try and incorporate the important dis- 
tinction between total void space and the available void 
space into the recent theories concerning compaction and 
packing 

In conclusion, we have studied theoretically density re- 
laxation towards steady state in granular compaction us- 
ing a microscopic model in one dimension. Due to vol- 
ume exclusion, exponentially growing time scales are as- 
sociated with cooperative motion of grains. As a result 
the approach towards the steady state is an inverse log- 
arithmic one. Since the argument leading to the loga- 
rithmic relaxation is a general one, the results should 
hold in a large class of physical situations, for example 
when the shaking is horizontal rather the vertical or when 
the grains are aspherical. Interestingly, the same mecha- 
nism is responsible for the long time necessary for packing 
grains in a bowl, parking a car in a busy street, or even 
entering a crowded room. 
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